Coupling Josephson qubits via a current-biased information bus 
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Josephson qubits without direct interaction can be effectively coupled by sequentially connecting 
them to an information bus: a current-biased large Josephson junction treated as an oscillator with 
adjustable frequency. The coupling between any qubit and the bus can be controlled by modulating 
the magnetic flux applied to that qubit. This tunable and selective coupling provides two-qubit 
entangled states for implementing elementary quantum logic operations, and for experimentally 
testing Bell's inequality. 

PACS. 74.50.-|-r - Proximity effects, weak links, tunneling phenomena, and Josephson effects. 

PACS. 03.67.Lx - Quantum computation. 

PACS. 03.65.Ud - Entanglement and quantum nonlocality (e.g. EPR paradox. Bell's inequalities, 
GHZ states, etc.). 



Superconducting circuits with Josephson junctions of- 
fer one of the most promising candidates for realiz- 
ing quantum computatio n [3, H, USHIESIEIS 
Efl 111 III III 111 111 HIT These superconducting 
qubits can be either charge- 01, flux- Q, mixed- H, 
current-biased Josephson-j unction (CBJJ) qubits Hlg, 
and others. Much attention is now devoted to realiz- 
ing controlled couplings between superconducting qubits 
and irnplementing quantum logic operations (see, e.g., 
011 ll in Il3)- Two qubits, i and j, can be con- 
nected by a common inductor or capacitor, with Ising- 

OTz). However, in 



a^J^ (a 



type couplings cr,„ 
general, (1) the capacitive coupling between qubits 

is not tunable (and thus adjusting the physical parame- 
ters for realizing two-qubit operation is not easy), and 
(2) a large inductance is required in to achieve a 
reasonably high interaction strength and speed for two- 
qubit operations ^1|- Alternatively, other schemes (see, 
e-g-, j^mOl) use sequential interactions of individ- 
ual qubits with an information bus. These provide some 
advantages: allow faster two-qubit operations, may pos- 
sess longer decoherence times, and are scalable. These 
schemes are similar to the techniques used for trapped 
ions |l7j |. where the ions are entangled by exciting and 
de-exciting quanta (data bus) of their shared collective 
vibrational modes. 

Compared to the externally-connected LC-resonator 
used in Rcf. ll and the cavity QED mode proposed in 
Ref. [U], a large (e.g., 10 ^im) CBJJ glU is more suit- 
able as an information bus, because its eigenfrequency 
can be controlled by adjusting the applied bias current. 
In fact, such data bus to couple distant qubits has been 
proposed in However, there all non- resonant inter- 
actions between the qubits and the bus were ignored. 
This is problematic because these near-resonance inter- 
actions must be considered, otherwise, the desired cou- 
pling/decoupling between the chosen qubit and the bus 
cannot be implemented because a perfect resonance con- 



dition is not always achievable. Also, modulating the 
bias current to selectively couple different qubits changes 
the physical characteristics (e.g., eigenfrequency) of the 
bus, and thus may yield additional errors during the com- 
munication between qubits. Finally, an effective method 
still lacks for refocusing the dynamical-phase shifts of the 
qubits to realize the desired quantum operations. 

Here, we propose an effective scheme for coupling any 
pair of superconducting qubits without direct interaction 
between them by letting these be sequentially connected 
to a large CBJJ that acts as a data bus. The qubit 
in Ref. [ij is a CBJJ, while here we consider charge 
qubits. Here, a large CBJJ acts only as the information 
bus between the qubits. Also, in contrast to Ref. [lTI |. 
in the present circuit any chosen qubit can be coupled 
to and decouple from the bus by switching on and off 
its Josephson energy. The bias current applied to the 
bus is fixed during the operations, and the dynamical- 
phase shifts of the qubits can be conveniently refocused 
by properly setting the free-evolution times of the bus. 
Therefore, an entanglement between distant qubits can 
be created in a controllable way for realizing quantum 
computation, and also for testing Bell's inequality. Its 
experimental realizability is also briefly discussed. 

Model. — Without loss of generality, we consider the 
simplest network sketched in Fig. 1. It can be easily 
modified to include arbitrary qubits. Each qubit con- 
sists of a gate electrode of capacitance Cg and a single- 
Cooper-pair box with two ultrasmall identical Joseph- 
son junctions of capacitance cj and Josephson energy 
ej, forming a superconducting quantum interference de- 
vice (SQUID) ring threaded by a flux $ and with a gate 
voltage V. The superconducting phase difference across 
the fcth qubit is represented by (jjk, k — 1,2. The large 
CBJJ has capacitance Cb, phase drop 0h, Josephson en- 
ergy Eh, and a bias current lb- The qubit is assumed to 
work in the charge regime with ksT <^ Ej <^ Ec ^ A, 
wherein quasi-particle tunnelling or excitations are effec- 
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FIG. 1; A pair of SQUID-based charge qubits, located on 
the left of the dashed line, coupled to a large CBJJ on the 
right, which acts as an information bus. The circuit is divided 
into two parts, the qubits and the bus. The dashed line only 
indicates a separation between these. The controllable gate 
voltage Vfc (k — 1, 2) and external flux $fc are used to manip- 
ulate the qubits and their interactions with the bus. The bus 
current remains fixed during the operations. 



tively suppressed. Here, fc^, T, A, Ec, and Ej are the 
Boltzmann constant, temperature, superconducting gap, 
charging and Josephson coupling energies of the qubit, 
respectively. The present mechanism of quantum ma- 
nipulation is significantly different from Refs. [^[lolllll|. 
although the circuits appear to be similar. The Hamilto- 
nian for the circuit in Fig. 1 is 
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with [4>k,nk] = i- Here, Qb — iirpb/^a is the opera- 
tor of charges on the CBJJ and [0f,,pb] 



2e[p cos(^$fc/$o),Cfc 



= Vfe/(2e), and Cb = Cb + ELi CyC^/Ck- 
$0 = h/(2e) and Uk are the flux quantum and the excess 
number of Cooper pairs in the superconducting box of the 
kth qubit, respectively. When the applied gate voltage 

Vk is set near the degeneracy points {ug''^ = (/ + 1/2), I = 
0, 1, 2, ...), then only the two lowest-energy charge states, 
\nk = 0) = I tfc) and \nk = 1) = | ifc), play a role. 
The large CBJJ works in the phase regime and describes 
the motion of a "particle" with mass m = Cf,(<i>o/27r)^ 
in a potential U{<f>b) = —Eb^coscjjb + It'f'b/Ic) with Ic — 
27r£'ft/<i>o. For lb < Ic, there exist a series of minima 
of U{<j)b)- Near these points, U{(f>b) approximates a har- 
monic oscillator potential with characteristic frequency 

ey{2Cb). 
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iUb = JsE^c Eb/f^^ [1 - iib/ic)^r\E: 
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/fc2 r-, IT ^2ll/4 

The approximate number of metastable quantum bound 
states Ha is TV, 



2^/^jEb/EP{l - h/Icf'^. For a 



restricted [yf to the Hilbert space spanned by the two 
lowest states of this data bus: |Ofc) and |lf,). 

The eigenenergy TvjJlc of the bus used in is much 
higher than that of the qubits. Therefore, adiabatically 
eliminating such a bus yields a direct interbit coupling. 
However, the energy scale of our proposed data bus (i.e., 
the CBJJ oscillator) is 0Jb/2TT ^ 10 GHz "e^l, which is of 
the same order of the Josephson energy (e.g., Ej/h ^ 13 
GHz 0). Therefore, the quanta in the present bus can 
be excited or de-excited when the qubit is operated. The 
Hamiltonian (1) clearly shows that the coupling between 
the chosen fcth qubit and the bus can be turned on and 
off 01, when the threaded flux $fc differs from or equals 
to (/' l/2)$o, I' = 0,1,2,.... For simplicity, hereafter 
we let Z, /' = 0. Two qubits can be indirectly coupled 
by independently interacting with the bus sequentially 
when exciting/de-exciting the vibrational quanta of the 
bus. Under the usual rotating-wave approximation, the 
dynamics for such a coupling mechanism can be described 
by the following effective Hamiltonian 



Hkb — Hk + Hb + iXk 



(2) 



Hk — 



with d = 
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■a^, Hb = hujb 
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and n — a^d being the Boson operators of the bus. Here, 



2e\l - 2nf>)/Ck, Xk - Cfc cos(7r$fc/$o), Cfc 



tors 



2h/{Ck'^oVCi 



b^b , 



The pseudospin opera- 



lfe)(lfe| - |0fe)(0fe|, = |lfc)(Ofc|, and 

CTl ' = |Ofe)(lfc| are defined in the subspace spanned by 
the logic states: |Ofe) = (life) -I- | Tfe))/\/2 and |lfc) = 
(life) — |Tfc))/\/2- Only the single-quantum transition 
process, approximated to first-order in ^f,, is considered 
during the expansion of the cosine-term of the Hamil- 
tonian (1), as the fluctuation of (j)b is very weak. In 
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fact, ci'^^i^D/Ck < 10 
imental parameters 0j S 



(k) /Mk) 



GHz, and Q VC, 



J 



10- 



<^ 1, for typical exper- 
Cb ~ 6pF,uJb/2Tr - 10 
Once the bias current lb 



is properly set up beforehand, various dynamical evolu- 
tions can be induced by selecting the applied flux $fe 
and the gate voltage Vk- Considering two extreme cases, 
the strongest coupling ($fe = 0) and the decoupling 
(<f>fe — $o/2), several typical realizable evolutions de- 
duced from the Hamiltonian (2) are given in table I. 



There, HAk = Sfe - fmjb, £k 



Jkh 



Quantum gates. — The physical characteristic (e.g., the 
eigenfrequency) of the bus in the present circuit does not 
need to be changed, once it is set up beforehand. It still 
undergoes a free evolution Uo{t) ruled by a non-zero Hb 
during the operational delay, i.e., the time interval when 
the qubits do not evolve because their Hamiltonians are 



low bias current, the dynamics of the CBJJ can be safely temporarily set to zero (when <I>fe = <&o/2, Vk = e/C, 



{k)^ 
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Controllable Parameters 


Evolutions 


Vk=e/Cf\ $fc = 'l>o/2 


t/o(t) 




ijf\t) 











TABLE I: Typical settings of the controllable experimen- 
tal parameters (Vfc and $fe) and the corresponding time 
evolutions tjj{€) of the qubit-bus system. Here, C'f'^ and 
le^j^ are the gate capacitance and the maximal Josephson 
energy of the feth SQUID-based charge qubit. (^^ is the 
maximum strength of the coupling between the fcth qubit 
with energy e-k and the bus of frequency lOh- The detuning 
between the qubit and the bus energies is /lAfc = Sh — fwjb- 
n — 0, 1 is occupation number for the number state |n) 
of the bus. The various time-evolution operators are: 
Uo{t) = exp{-itHi/h),U[''\t) = exp[-it(5S^'='ai'''/(2ft)l ® 
Uo{t),u!t^ = i(t)cosA„|0fc)(0fc|-(sinA„)a/V7rrTl0fe)(lfcl + 
atsin|„/y^|U)(Ofc| + cos|„!0fc)(0fcl, and U^^^t) = 
i(t)exp{-itCi[|lfc)(lfc|(ri + 1) - |Ofc}(Ofc|w]/(fiAfc )}, w ith 
A{t) = exp[-if(2Ji'6 + E''j^a^J'^)/{2h)], A„ = 2CktyM+l/h, 
and 5„ = 2C,ktVfi/fi. 



Before and after the desired operations, the bus should 
remain in its ground state |0f,). In principle, the time- 
evolutions listed in table I are sufhcient to implement 
any desired operation for manipulating the quantum in- 
formation stored in the present circuit. In fact, any 
single-qubit rotation, including the typical Hadamard 



gate: H, 



(k) _ 



]/V2, on the fcth qubit can 



be easily realized by selectively using ?7f)(t) andJ/f^i)- 
Any single-qubit operation is not influenced by the free 
evolution of the bus during the time delay between suc- 
cessive operations. 

In order to realize two-qubit gates, we must be able to 
couple distant qubits via their sequential interactions to 
the bus. We set the bias current h such that hwh = 2e^P 
and then perform a three-step process. First, we cou- 
ple the first (control) qubit to the bus by switching off 
its appHed flux <I>i and produce the evolution {j^\ti). 
After a duration ti determined by sn\{2(^iti/K) — —1, 
the control qubit is decoupled from the bus. This pro- 
cess implements the evolutions: |0f,)|0i) |0&)|0i) and 
|Ofe)|li) e-'"''*i|l6)|0i). Next, we let the second (tar- 

f 2") 

get) qubit work at a non-degenerate point (V2 7^ e/Cg ') 
and couple it to the bus by switching off its applied flux 
$2- After a duration ^2 determined by 



sm 



(i 



2I2 



cos 



£2^2 

2h 



C¥2 



= 1, 



(3) 



the target qubit is backed to its degenerate point (V2 = 

f2) 

e/Cg ) and decoupled from the bus. This sequence of 



operations generate the evolutions: 

|0t)|02) ^ e-'«|0fc)|02), |0b)|l2) e-^«|0b)|l2), 
I Ifc) IO2) ^ i e-''(?+-.t.) (cos 7^2 1 U) IO2) + sin 772 1 U) 1 12)) , 
|lh)|l2) ^ie-^(«+"''*»Usinr72|lfc)|02)-cosr?2|U)|l2)), 

with ^ = ujb{Ti +t2 + T2)/2 + (it2/{2h^^2), and tg = 
X]s=i ^ X]s=i '^s- Finally, we couple again the control 
qubit to the bus and perform the evolution {t^) with 
sin(2Cii3//i) = 1, yielding evolutions: |06)|0i) |0h)|0i) 



and |lb)|Oi) 



|Ofc)|li). In practice, the free evo- 



lutions Uq{ti) and Uq{t2) exist during the time delays be- 
tween the first (second) and second (third) pulses. If the 
delays are further set accurately such that the total dura- 
tion tg satisfies the condition sin ujbtg = 1 , then the above 
three-step process with two delays yields a quantum op- 



eration U{tg) = u!2\h)Uo{T2)Ul^> {t2)Uo{n)U^'> {t 
exp(~iO|0,,)(06| ® Uf^\(32). with 

/l 



(2), 



\ 

10 

cos/32 sin/?2 

sin /32 — cos /32 / 



(4) 



(2'] 

being a universal two-qubit gate. Here, cos/32 = 2ej /e2- 
This gate can produce entanglement between qubits and 
also realize any quantum computation, accompanied by 
single-qubit rotations. 

Testing Bell's inequality. — Entanglement is a key in- 
gredient for computational speedup in quantum compu- 
tation. Historically, Bell's inequalities were seen as an 
entanglement test: its violation implies that entangle- 
ment must exist. For a two-qubit entangled state |'i/'e), 
the Clauser, Horne, Shimony and Holt (CHSH) form of 
Bell's inequality : /(IV'e)) < 2 is usually tested by ex- 
perimentally measuring the CHSH function /(|'0e)) = 
|£;(0i , 02) + E{e[ ,62) + E{9i ,e'^)-E{e[,9'2)\. Here, 9^ are 
controllable classical variables and E{6i, O2) is the corre- 
lation for the outcomes of separately projected measure- 
ments of two qubits. A number of experimental tests poj 
of Bell's inequality have already been performed by using 
entangled photons and atoms. We now show that a de- 
sired entangled state can be created in a repeatable way 
and thus Bell's inequality can also be tested experimen- 
tally by using this circuit. 

We begin with an initial state — \Qb)\ii)\i2!^\Qb)® 
(|0i) + |li)) (g) (IO2 + II2)) /2 with two qubits decoupled 
from the bus but working at their non-degenerate points 

(i.e., Vk ^ e/Cg'^^). After applying a Hadamard gate 

^ (2) 

Hg ' to the second qubit, the system evolves to the state 
= |Ofc) ® (|0i) + 111)) ® |l2)/V2. The desired two- 
qubit entangled state is then generated as 



P2))^u['\9,)Uri92)Urm\^i), (5) 



r(2 



(12), 



with Ul'''i9k) = exp i9kai'''^/2 



5E^^hk/h. The 



corresponding correlation function is E(9i,92, (32) 
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sin (sin 6*1 sin 6*2 — sin /32 cos cos 62 ) ■ For certain 
chosen sets of angles: 9k — {— 7r/8, Svr/S}, the CHSH 
function becomes 



\/2 I sin/32(sin/?2 



(6) 



It is easy to numerically check that Bell's inequality, 

f{\^T\P'i))) < 2, is violated when iJ^^V^^c^ < 0.776, 
which can be easily satisfied for this charge-qubit system. 

Experimentally, the above procedure can be re- 
peated many times at each of the four sets of an- 
gles and thus the correlation function E(9i,02, P2) — 

[^same(^l,e2) - A^diff (^1 , ^2)] Mot , with iV^amc (^1 , ^2 ) 

(.^diff(^i, ^2)) being the number of events with two qubits 
being found in the same (different) logic states and 
A^tot = A^samc(6'i,6'2) + iVdiff (6*1 , 6*2) being the total ex- 
perimental times for the same 9i and 02. Finally, Bell's 
inequality can be tested by calculating the experimen- 

tal CHSH function: f{\^T\f32))) = \E{ei,e2, P2) + 

Ei9[,02.P2) + Eiei,e'2,P2)-E{e'ye'2,i32)\. 

Discussion. — Two types of noise, fluctuations of the 
applied gate voltage 14 and bias current h, must be con- 
sidered in the present qubit-bus system. For simplic- 
ity, these two environmental noises are treated as two 
separate Boson baths with Ohmic spectral densities and 
assumed to be weakly coupled to the qubit and CBJJ, 
respectively. The Hamiltonian of the fcth qubit coupling 
to the bus, containing these fluctuations, can be written 
as 



eC, 



(fc) 



H = Hkt+ E 

3 = 1,2 ^ 



2CbUJb 



a^R2 + aRl \ , Rj gu^a^, ■ (7) 



decoherence rates are characterized [2l| by two dimen- 
sionless couphng parameters, ay = {Cg'^^ /Ck)'^ Rv / Rk 
and ai = Re(Yj)/ {Cb^^b), which describe the couplings of 
the voltage fluctuations to the qubit and the bias-current 
fluctuations to the bus, respectively. Here, Rk = h/e"^ 
25.8 kfi is the quantum of resistance, Rv is the Ohmic 
resistor of the voltage, and Re{Yj) is the dissipative part 
of the admittance of the current bias. If the qubit decou- 
ples from the bus, ay (ai) characterizes the decoherence 
and relaxation of the qubit (bus). It has been estimated 
in Ref. Jj that the dissipation for a single SQUID-based 
charge qubit is sufficiently weak (ay « 10^®), which al- 
lows, in principle, for 10^ coherent single-qubit manip- 
ulations. However, for a single CBJJ the dimensionless 
parameter aj only reaches 10^'^ for typical experimen- 
tal parameters 1^: l/Re(F/) ~ 100 n, Cb ~ 6 pF, 
bJb/^n ~ 10 GHz. This implies that the quantum co- 
herence of the present qubit-bus system is mainly lim- 
ited by the bias-current fluctuations. Fortunately, the 
impedance of the above CBJJ can be engineered to be 
l/Re(Yf) - 560 Ml "ej. This lets a/ reach up to 10~^ and 
allow about 10^ coherent manipulations of the qubit-bus 
system. 

In summary, we have proposed a scheme for coupling 
two SQUID-based charge qubits by sequentially using 
their interactions with a common large Josephson junc- 
tion biased by a fixed current. Each interaction is tun- 
able by controlling the external flux applied to the cho- 
sen SQUID-based charge qubit. The proposed circuit al- 
lows the possibility of implementing elementary quantum 
logic operations, including arbitrary single-qubit gates 
and universal two-qubit gates. The created two-qubit 
entangled states can be used to test Bell's inequality. 



Here, Qui^ is the Boson operator of the jth bath and Quj^ its 
coupling strength. The relaxation and decoherence rates 
of our qubit-bus system can also be calculated by using 
the well established Bloch-Redfield formalism [l^. Un- 
der the usual secular approximation, the relaxation and 
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